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Abstract 

The radiative decay of the massive neutrino i>i — > Vjj is investigated in the framework 
of the Standard Model in external electromagnetic fields of various configurations: con- 
stant crossed field, constant uniform magnetic field, plane monochromatic wave's field. 
The effect of significant enhancement of the neutrino decay probability by the external 
field (electromagnetic catalysis) is discussed. An especially strong enhancement occurs 
in the case of the ultrarelativistic neutrino decay, since in this case the decay probability 
does not contain suppression caused by the smallness of the decaying neutrino mass. The 
ultrarelativistic neutrino decay catalysis is significant even in a relatively weak external 
field (F/F e <C 1, where F e is the critical Schwinger value). The expression for the photon 
splitting probability into the neutrino pair 7 — * ViVj in the wave field is given. The estima- 
tions of a number of gamma-quanta produced in a volume filled with an electromagnetic 
field and the neutrino lifetime in a strong magnetic field are presented. 
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1 Introduction 



For quite a long time we have seen intensive theoretical studies of flavour- 
changing processes caused by the phenomenon of fermion mixing. The descrip- 
tion of this phenomenon in the quark sector goes back to the pioneer work by 
Cabibbo |I| and presently is put into practice by introducing the unitary 3x3 
matrix Vfj (the so called Cabibbo-Kobayashi-Maskawa matrix ||. It should be 
noted that qualitative progress has been observed in experimental investiga- 
tion of quark mixing parameters. Except a detailed examination of the "old" 
mixing angles related to the first two quark generations, important information 
has been obtained from the decays of 6-quark containing particles (ARGUS 
Collaboration |3j], CLEO Collaboration |§]). On the other hand, thus far there 
is no experimental evidence in favor of the analogous mixing phenomenon in 
the lepton sector. This can be accounted for in a natural way by the fact that, 
because of the insufficiently high precision achieved in experimental studies of 
neutrino-involving processes, the neutrino mass spectrum appears degenerate 
(the neutrinos manifest themselves as massless particles j5j]). The neutrino 
mass spectrum being degenerate, lepton mixing is known to be purely formal 
and unobservable. At the same time, with massive neutrino the absence of 
lepton mixing seems unnatural and is virtually incompatible with attempts to 
somehow extend the standard model. Notice that lepton mixing may lead to 
some interesting physical phenomena such as: 

1) charged lepton radiative decays with lepton number violation of type fi — ► ej, 
/i^3ei,0,/i^e 7 7[i. 

2) neutrino radiative decays v\ — > ||, V{ — ► z^77 0; 

3) neutrino oscillations [|H| , 

4) the possible effect of massive neutrino mixing on the spectrum of (3 - decay- 
produced electrons [|TTj. 

Even such a short review of lepton mixing effects shows that most of these 
are associated with the massive neutrino. Nowadays, the physics of the massive 
neutrino is becoming a vigorously growing and prospective line of investigation 
at the junction of elementary particles physics, astrophysics and cosmology. It 
will suffice to mention the well known problem of the solar neutrino [jT2j] and the 
possibility of solving it (the mechanism of resonance enhancement of neutrino 
oscillations in substance the effect of the massive neutrino radiative decay 
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on the spectrum of the relic radiation |T4| and so on. The above-mentioned way 



of solving the solar neutrino problem using the MSW mechanism^ shows that 
the the massive neutrino's properties are sensitive to the medium it propagates 
through. Substance is usually considered as medium. We note, however, that 
medium can also be represented by an external electromagnetic field, which 
can significantly influence both the properties of the massive neutrino itself §T6 



and the process of its decay |T7j and even induce novel lepton transitions with 



flavour violation v\ <-> Vj (i ^ j) (jTSj], forbidden in vacuum. In our preliminary 
communication (jT7j we have pointed out the probability of the massive neutrino 
radiative decay V{ —> i/jj (i j) being considerably enhanced in a constant uni- 
form magnetic field. Such an enhancing influence of an external field can be 
illustrated with the straightforward example of a neutrino radiative decay in 
a weak (as compared with the Schwinger value F e = m\je ~ 4.41 ■ 10 13 G) 
electromagnetic field. To this end we use the amplitude of the Compton-like 
process Uff* — ► za,7* with virtual photons which, in particular, allows 



obtaining the first term of the expansion of the radiative decay V{ — > vf^ ampli- 
tude in a weak external field. In the expression for the amplitude of the process 
z/ ?Cpi)+7*( ( /i) ~~ *■ z/ j(j 9 2)+7*(^2) it is sufficient to consider 7(92) as a real photon, 
and to replace the field tensor of the virtual photon 7* (51) by the Fourier image 
of the external electromagnetic field tensor. Below we shall give the expression 
obtained in this way for the radiative decay amplitude in the simplest case of 
a uniform electromagnetic field, in which the decay kinematics p\ + = P2 + #2 
is the same as in vacuum. The external-field-induced contribution AA4 to the 
amplitude of the decay — > ^7 can be represented in the following form: 

where = ^(^2)7^(1 +75)^(^1), i,j = 1,2,3 enumerate the definite mass 
neutrino species, pi, P2, Q are the four-momenta of the initial and final neutrinos 
and the photon, respectively, is the external uniform electromagnetic field 
tensor, f a p{q) = ^a/3^q^v{q), £v{o) is the polarization four- vector of the photon, 
Ft = mj je is the critical value of the strength of the electromagnetic field for the 
charged lepton with the mass mi. We have introduced the following designation 



1 About the current situation around the solar neutrino problem see, for example, ref's [15 
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< A(m e ) >= £ K ie K* e A(m e ), (2) 

£=e,[i,T 

where Km (£ = e, /i, r) is the lepton mixing matrix of Kobayashi-Maskawa type. 
For the sake of comparison we write the known expression for the amplitude of 
the neutrino radiative decay in vacuum M which can be represented as follows: 



where p = pi + P2, me, mw are the masses of the virtual lepton and VF-boson, 
respectively. In analyzing the amplitudes ([I]) and (0) in the case of the neutrino 
decay at rest, it is necessary to take account of p%, p^-, q, j being of order of 
the mass of the decaying neutrino m v . In this case, the expressions for the 
amplitudes ([]]), (§) can be easily estimated (it is sufficient to allow for the order 
of the dimensional quantities): 



AM ~ G F ml(F/F e ), (4) 
Mo ~ G F ml(m T /m w ) 2 . (5) 

It follows here from that, given the condition 

(F/F e ) 2 » (m T /m w ) A (6) 

(here F stands for the strengths of the magnetic (B) and electric (E) fields), the 
probability of the decay V{ — > vfj in an external field is much greater than that 
in vacuum, even for a relatively weak electromagnetic field (10~ 3 <C F/F e <C 1). 
The catalyzing effect of an external field becomes even more substantial in the 
case of the ultrarelativistic neutrino decay [E v ^> m v ). With the amplitudes 
(S)? (0)? being Lorentz-invariant, the analysis can be conveniently carried out 
in the rest frame of the decaying neutrino. In this case the the electromagnetic 
field in Eq. (£|) is obtained by the Lorentz transformation from the laboratory 
frame, in which the external field F is given, to the rest frame of the decaying 
neutrino: 

F' - —F > F. (7) 

TTly 
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Comparing the expressions (g) and (|5j), in view of (0), we notice that the cat- 
alyzing effect of the external field becomes appreciable under a much weaker 
condition, as compared to (§): 



(pxFFpt) (m T \ 



) ■ (8) 



raj; F e 2 \mw, 

In this case the ratio between the probabilities of the ultrarelativistic neu- 
trino decay and the decay in vacuum is of the order: 



The expression (|9|) shows that in the ultrarelativistic neutrino decay the en- 
hancement is mainly due to a decrease in the decay probability suppression 
by the smallness of the neutrino mass (w^ ~ m\, ~ m h v (m v / E v )) . It is 
natural to expect that in taking the account of further terms in the expansion 
of the amplitude of the radiative decay V{ — ► vfi with respect to the external 
field, the suppression mentioned above can be fully canceled. All this makes 
it interesting to calculate the amplitude with the external electromagnetic field 
taken into account exactly. An expression thus obtained will be valid in the 
case of the neutrino radiative decay V{ — ► vfj in an external electromagnetic 
field, which has not to be weak as against the Schwinger value F e . 



2 The crossed field 

At present the experimentally accessible strengths of electromagnetic fields are 
significantly below the critical strength (F/F e <C 1, F = B,S, F e = m 2 e /e ~ 
4.41 ■ 10 13 G). Because of this , field-induced effects are especially marked in 
the ultrarelativistic case with the dynamic parameter 

x 2 = ^ (10) 

777° 

being not small even for a relatively weak field (F^ is the external field tensor, 
p a is the 4-momentum, m is the mass of the particle). This is due to the fact that 
in the relativistic particle rest frame the field may turn out of order of the critical 
one or even higher, appearing very close to the constant crossed field. Thus, 
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— * — * 

the calculation in constant crossed field [S _L B, 8 = B) is relativistic limit of 
the calculation in an arbitrary weak smooth field, possesses a great extent of 
generality and acquires interest by itself. We note that, as (FF = FF = 0) 
in crossed field, the dynamic parameter x 2 (033) i s the single field invariant, 
by which the decay probability is expressed. Furthermore, the calculation in 
crossed field is the least cumbersome and, therefore, we consider this case first 
just to outline the calculation technique. 

In the lowest order of the perturbation theory, a matrix element of the ra- 
diative decay of the massive neutrino V{ — ► vff (i ^ j) in the Feynman gauge is 
described by the diagrams, represented in Fig.l (a,b) where double lines imply 
the influence of the external field in the propagators of intermediate particles. 
Summation is made over the virtual lepton £ in the loop (£ = e,/i, r). Un- 
der the conditions mf/m^r <C 1, eF/rriyy <C 1 the field induced contribution 
AA4^ = M. — to the decay amplitude can be calculated in the local 
limit, in which the lines W and ip are contacted to a point. It is most easily 
seen, if &MS F > is expanded into a series in terms of the external field, as is 
shown in Fig 2, where the dotted lines designate the external electromagnetic 
field A ex . We note that the first seven diagrams in Fig 2 coincide with the 
diagrams describing Compton-like process vf<j* — > v^* and, as was pointed out 
nUl , this process amplitude is reduced to the contribution of the first two 



m 



diagrams in the local limit. With the orthogonality of the mixing matrix 
taken into account, this is due to the fact, that the main contribution to the 
integral over momentum in the loop gives from the region of the virtual mo- 
menta p ~ mi <C raw- We remind that we are investigating flavour violating 
processes (i ^ j), and, hence, < A >= 0, if A is independent of mi . Thus, 
the dominant contribution of order l/m^ ~ Gp only comes from the diagrams 
with one VF-propagator in Fig. 2. 

Even such a simple analysis shows the following: 

1) the predominant contribution to the amplitude is made by the diagram (a) in 
Fig.l. This diagram in the local limit of VF-boson propagator being contacted 
to a point, transforms to the diagram shown in Fig. 3; 

2) since in calculating A.M the mass of the VF-boson in the local limit appears 
in the weak interaction constant Gp = g 2 /8m^ v only, the amplitude does not 
contain the known GIM suppression factor of the decay U{ — ► Vf] in vacuum 
- m\lm 2 w <C 1 (see Eq.(g)). 
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The expression for the amplitude, corresponding to the diagram in Fig. 3, 
can be represented in the following form: 



AM = 



ieGp 

-j=- 3i<-<A<l) <* ■'<>■ i 



e* a {q) < J a p(q) > -M (0) , 



d 4 XSp j a S(X) lf3 (l + l5 )S(-X) 



12) 



where X = x — y and S(X) is the propagator of a charged lepton in the crossed 
field (see Appendix A, Eqs.(A.l) and (A.2)). All the other quantities in ( |TT| ) are 
defined above (see Eqs.(|Tj)-(@)). The details of the tensor J a p{q) calculation may 
be found in Appendix A, while here we only give the result of the calculation: 



AM 



(qFFq) Smj 



+ J^{Fh{qF3)h + -^-AFr){qF3)h), 



24mj 

fa (3 = qa£(3 — qpZ 
1 



48m! 



(13) 



ai 



'\WOL 



faf3i 



where F^, F^ u = e^ a p F a p/2 are the tensor and the dual tensor of a constant 
field; e > is the elementary charge, Gp is the Fermi constant. In Eq . (p~3|) I a 
(a = 1, ... ,4) are integrals of the known Hardy-Stokes functions f(u): 



h 

h 

h 

h 

m 



dt uf(u), 
j l Q dt{l-t 2 )uf{u), 

f'dtil-t 2 ) (3-t 2 )u 2 ^-, 
Jo v ' v ' du 

f'dtil-t 2 ) (3 + *V^, 
v 1 v 1 du 

1 

dz exp 



% (zu H — z 3 ) 
v 3 7 



(14) 



(15) 



u 



[e 2 (qFFq)/16ml)- 1 / 3 (l-t 2 )- 2 / 3 



As can be readily checked, the amplitude ( |T3| ) is evidently gauge invariant, as 
it is expressed in terms of the tensors of the external field F^ v and the photon 
field f^ v . In the weak field limit (F/F e <C 1) the predominant contribution to 
the amplitude is made by the second term in the epression (|T31 ). The integrals 
I a can be easily evaluated in this limit, taking into account the orthogonality 
of the lepton mixing matrix Ku : 



<h> = < 1 + 0[(F/F e ) 2 } > ~ 0[(F/F e )\ 

h ^ 2/3, (16) 

h - 28/15, 

h ~ 0[(F/F e )% 

As expected, the amplitude (|T3| ) in view of flI5|) in the weak field limit coincides 
with the expression ([]]). The amplitude of the process Vi — > in a crossed field 
([131) is simplified substantially in two cases: that of the decay of a neutrino at 
rest [E v = m v ) and that of the decay of an ultrarelativistic neutrino [E v ^> m v ). 



2.1 Neutrino at rest {E v = m v ) 

In this case, the dynamic parameter (|ID|) 

2 _ e 2 (piFFpi) _ lm v eF\ 

— 6 — 2 

rrq \ rri£ mf J 

is obviously small even when the field strength exceeds the critical values (F > 
mj/e, m v <C rri£, %t ^ !)• The decay amplitude ( |I3|) , (|TJ) can then be reduced 
to the form: 



~ ^n(qFj)}(KieK* e ). (17) 

Here we have introduced the dimensionless field tensor J 7 ^ = F^/F e: where 
F e = m? e /e is the critical strength value, m e is the electron mass. It is clear 
from Eq.(|T71) that the decay probability is represented by a polynomial of the 
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sixth degree in field strength. In the limit F « F e (a weak field limit), the 
expression for the decay probability is determined by the lowest power of F 
and has the following form: 



„. t ~ _£l -9L m ? I x _ "AY (LX \ Ki k> p (is) 

~ 18tt 1927T 3 * \ ml) \FJ ' J e ' ' 



where and mj are the masses of the initial and final neutrinos. In the 
opposite case, F ^> F e (a strong field limit), we have: 

a Gl , / m5\ 5 / m 2 A /F\ 6 , „ l9 



The expressions and (J9j) should be compared with the well known prob- 




ability of the decay V{ — > z/?7 in vacuum 

1 A m > ( 1 + ^J l 1 "^ W (20) 

The comparison demonstrates the catalyzing effect of the external crossed field 
on the decay probabilities, as there is no suppression ~ (mf/mj^) in Eqs.([T8|) 
and Actually, the enhancement influence of the external field takes place 
in the weak field limit (F <C F e ) under the condition: F > 2 ■ 10~ 3 F e . Besides, 
the decay in the strong crossed field ( [191) is catalyzed by an additional factor of 
the form - (F/F e ) 6 > 1. 



2.2 Ultrarelativistic neutrino {F v ^> m u ) 

Notice that in the ultrarelativistic limit the kinematics of the decay Vi(pi) —> 
Vj{p2) + 7(<?) is such that the momentum 4- vectors of the initial neutrino p\ and 
the decay products p2 and q are almost parallel to each other. Therefore, the 
current 4-vector j a = ^ (£>2 ) 7a (1 + 75)^(^1) is a l so proportional to these vectors 
(ja ~ Pia ~ Qa ~ P2a) ■ In this case, the expression for amplitude (]T3j) can be 
simplified and reduced to the form 

1 

2 </i>, (21) 



AM 



e 2 G 1 



7T- 



e*F Pl 



[l-x) + 



l+x) 
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where x = cost?, d is the angle between the vectors p\ (the momentum of 
the decaying ultrarelativistic neutrino) and q 1 (the photon momentum in the 
decaying neutrino z/j rest frame). The argument u of the Hardy-Stokes function 
f(u) in the integral I\ (see Eq.(|PD) in the ultrarelativistic case has the form: 



u = 4 



l + x)(l-t 2 ) 1 



; 2 2) 



The Lorentz-invariant decay probability W-Ej, can be expressed in terms of 
the integral of the squared amplitude over the variable x: 



wE„ ~ 



167T 



m- 



dx \AM\' 



(23) 



a G 2 F 6 2 

^ -i m eX e 




'l-aO + 



1 + 



</i> 



For small values of the dynamic parameter (xi ^ 1) 5 the integral h(xt) is 
expanded into the following series: 



h 

Xi 



1 + + 4^ + 

1 + x ( . m 



j 



TTIa 



and the probability ( P3|) can be represented in the form: 



W ~ 



m 



4tt (15tt) 3 £ 

As the dynamic parameter ~ 
trino's energy, it is clear from Eq. 



1 



TTIa 



TTIa 



1 + 5- 



TTIa 



\K iP Kt\ 2 



(24) 



(25) 



[E v /mi){F /Ft) is proportional to the neu- 
f|25Q that, with increasing the energy of the 
decaying neutrino, the decay probability increases as ~ E h v . 

For great values of \i ^ 1? using the asymptotic behavior of the Hardy- 
Stokes function both at large and at small values of the argument and also the 
unitarity of the mixing matrix one can represent Eq. (|23|) in the form: 
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w _ G|m| 2 [ \KieK*f, Xe » 1, X,,r « 1, ^ 

47T 7T 3 ^ Ae \ mf J I |i^ r | 2 , Xe » X M » 1, Xr « 1, 

In this way, the decay probability increases linearly ~ E v (x e > 1, Xt ^ 1) 
and finally (x T ^> 1) becomes a constant: 

21.7a GJ, m? , r ^ |2 _ 

^ f -TT Xr (27) 

(TV?- \ 
1 " -| > 1. 



Comparing the decay probabilities ( |25|) -(|27|) in the crossed field with the vac- 



uum decay probability ( |2D|) , we notice that the catalyzing effect of the field on 
the ultrarelativistic neutrino decay becomes even more substantial compared to 
the situation with the neutrino at rest, because there is no suppression caused by 
the smallness of the mass of the decaying neutrino. Recall also that none of the 
epressions for the decay probability in the crossed field contain the well-known 
suppression GIM-factor (mi/rriw) 4 characteristic of the probability of the de- 
cay Vi — > in vacuum. The probability of neutrino decay at rest in a strong 
crossed field (see Eq . QTgj) ) is enhanced by the additional factor ~ (F/F e ) 6 ^> 1. 

Here we estimate a number of gamma-quanta which can be result as the 
decay product from a neutrino beam of a high energy accelerator in a volume 
filled with an electromagnetic field. As the experimentally accessible strengths 
of electromagnetic fields are significantly below the critical strength (F <C F e ) 
in the laboratory conditions we can use the expression (|Z5|) for the decay prob- 
ability of a high energy neutrino [E v ^> m e F/F e , Xe ^> !)• 

The number of gamma-quanta AN 1 which are produced in a volume filled 
with a magnetic field of the strength B transversal to the neutrino beam can 
be presented in the following form: 

^~ 10T S)\(^)(i5^)W' (28) 

where L is a longitudinal dimension of the volume. The parameter W has a 
simple meaning of a full energy of neutrinos passed through a "target" during 
the experiment. It can be presented as 
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where dN v ' jdE is the energy distribution of the neutrino beam. One can esti- 
mate this parameter from a number N cc of v' s charge current (CC) events in 
detector with known design parameters during the time of neutrino experiment. 

With the expected data on N cc from CERN-SPS neutrino beam which is 
presented, for example in ref. |2Uj], one can estimate: 



W ~ 10 19 GeF. 

As one can see from fl251), AN 1 ~ 1 requires B ~ 10 10 G in a volume ~ 1 m £ 



3 Constant Magnetic Field 

The probability of the massive neutrino decay V{ — > in a constant magnetic 
field having the strength B is described by two invariant parameters: 

1) the above-mentioned dynamical parameter 

2 _ e 2 ( Pl FF Pl ) _ B 2 p? 1 
X > m\ B|^| 1 (29) 

where Bi = mf/e is the critical magnitude of the magnetic field, pi± is the 
initial neutrino's momentum component, normal to the magnetic field; 

2) the external magnetic field intensity parameter 

, _ _e^(FF 1 _ & 

where F^ v is the constant uniform magnetic field tensor, (FF) = —2B 2 . Since 
magnetic fields in up-to-date superconductive magnets range up to strengths 
B < 10 5 G, the hope to obtain marked quantum effects induced by a magnetic 
field, as it seems, should be related only to the ultrarelativistic neutrino decay, 
when the dynamic parameter xe may be not small, while the intensity parameter 
r] e <C 1. In this case we have the crossed field limit which was considered in detail 
in the preceding section. The decay in a strong magnetic field (B > B e = 4.41 • 
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10 13 G) is probably of interest in astrophysics or in early Universe cosmology. By 
mentioning astrophysics, we, first of all, mean intensive magnetic fields, "frozen 
in" neutron star substance. The primordial magnetic fields, "frozen in" the 
cosmological plasma could also have had the strenghts B > B e at some stages 
of the early Universe evolution [j2T[ . 



The amplitude of the decay V{ —> Ujj in a magnetic field is described by the 
same effective Feinman diagram as that shown in Fig. 3. The expression for the 
lepton propagator and the calculation details are given in Appendix A, here 
we only present the external-field-induced contribution AA4^ to the decay 
amplitude. 

+ cfaXgw) < ^ > +i(fm^j) < n > ]■ (3i) 

We recall that here j M = 1^(^2)7^(1 + 75)^(^1) is the neutrino current, 

<Y{ mi )>= £ K u K* jt Y(mi), (32) 

Pi, P2-, q are the 4-momenta of the initial and final neutrino and the photon, 
respectively, f a p = q a €p — qpe a is the photon field tensor. We have introduced 
the dimensionless tensors of the external magnetic field cp a p = F a p/B, {p a p = 
tapLLvVuv/Z, A a/3 = ((p(p) a p- The double integrals Yi (i = 1, 2, 3, 4), entering 
(|3Tp , have the form: 



1 00 7 

az 



sin z 

1 — cos z cos zt 





l/i = (1 — t ) sinz + t sin zt 



sin z 



,. ?N tsmzt , ooN 

1)2 = cos zt — (1 — t ) cosz , (33) 

tan z 

t sin zt cos zt — cos z 

?/3 = — cos zt+- \-2 — * , 

tan z sin z 

z/4 = 2j/i - (1 - £ 2 )sinz, 



13 



n (2)t ) = - f£2-g) + J_-^ (£M_£. 

2 y 2 tan z smz y e5 77^ 

It should be noted that in ( g3p the integration with respect to 2 is performed 
over the complex 2; plane. The integrand's singularities are bypassed in a usual 
way |3Uj in the lower semi-plane (Im z < 0) of the complex z plane. In the weak 
field approximation B <C B^ {r\i <C 1) the integrands in fl33Q become rapidly os- 
cillating functions of the variable z (since at z ~ 1 we have £) ~ l/r]g^> 1). 
Therefore, the predominant contribution to the integrals comes from a narrow 
region of small z ~ 77^, in which the integrals can be considerably simplified. 
By straightforward calculation, it can be shown that the amplitude (pTj) in this 
limit (77^ = <C 1, however = (B/B£)(p±/mt) has not to be small), 

as should be expected, coincides with the amplitude of the neutrino radiative 
decay in the crossed field (|T3|). Note, that this is the necessary (but, certainly, 
not sufficient) condition for our combersome calculations of the magnetic-field- 
induced amplitude AA4^ to be correct. 

In the strong field approximation B ^> Bi (77^ ^> 1), to evaluate the inte- 
grals (|33|), it is convenient to rotate the integration loop clockwise by tt/2 in the 
complex z-plane. In this case the amplitude £±MS B > (PT|), (|53f) is significantly 
simplified and can be represented as follows: 



AMi? ^ ^ % (ffilUQ) + (m) + UM)} < mH^mVql >, (34) 

where q\ = (qAq) is the square of the photon momentum, which is normal to 
the magnetic field, 



H(x) = 



H(x) = 



3 

2 X 



x 



arctan 



x 



x 



X 



X 



In 



x 



v 

i + yr 
1 - vr 



1 



X > 1. 



(35) 



X 



X 



+ 2 + 17T 



X 



X < 1. 



X 



Note that, if q\ = 4mf (x = 1), the amplitude (fJ|) has a root singularity 
associated with the known root singularity of the probability of the e + e~ pair 
generation by a photon in an external magnetic field |j22fl . Such a behavior 
of the decay amplitude in the vicinity of q\ — > 4mf exactly corresponds to 
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the singularity of the imaginary part of the photon polarization operator in a 
magnetic field j23( . 

On the basis of the expression (^) for the amplitude we obtained the follow- 
ing expression for the probability of the ultrarelativistic neutrino of a moderate- 
energy {jp\ < 4m|) in a strong magnetic field (rjt ^> 1): 



2a Gl m 6 p ( B \ 2 , r , _ l2 T/ . 
J(z) = | dy(z-y) 



1 -i 2 





2 = 



1 y _ 

^—^arctan-^-l 



sin 6 



2m 



e 



To estimate the enhancement of the decay by the external magnetic field numer- 
ically, it is sufficient to compare the expression ( |35Q for the decay probability 
in a magnetic field with the expression flZUQ for the probability of the decay of 
the massive neutrino in vacuum. 

As we can see, in the relativistic neutrino decay probability ( |3B|) there is no 
suppression caused by the smallness of the neutrino mass, because the proba- 
bility is virtually independent of the neutrino mass (if m 2 /m 2 <C 1). Finally, 
the decay in a strong field is also catalyzed by a factor ~ (B/B e ) 2 . Especially 
impressive is the comparison of the moderate-energy [E 2 < Am 2 ) relativistic 
neutrino's lifetime 

T {B) - m 

the radiative decay V{ — > vfy in a strong magnetic 

field {B > B e ) 
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7 /» \ 2 



and the lifetime 

r (0) in 

the radiative decay in vacuum 
((]) 10 50 fleV\ 6 ( E v \ 

So, the magnetic catalysis of the massive neutrino radiative decay might solve 
the problem of whether a heavy enough {m v > 20 eV) neutrino exists in the 
Universe. In fact, if the magnetic field fluctuations in an early enough (t > 1 sec) 
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Universe would, for some reasons, reach such values as B ~ lO 3 ^, the massive 
neutrino's lifetime under such fields could, according to ([T7]), reduce to values 
of order of one second. 



4 Plane monochromatic wave 

The development of intensive electromagnetic field generation techniques and 
the current possibility to obtain waves of high strength of electromagnetic field, 
namely E ~ 10 9 V/cm, stimulate the investigation of quantum processes in 
strong external wave fields. Indeed, the so called wave intensity parameter 



2 / or\ \ 

&e = r> ( 39 ) 



171} 

(where a is the amplitude of wave, mi is the lepton mass and e is an elementary 
charge), characterizing the effect of the electromagnetic wave should not be 
neglected. It's also worth noting that the energy-momentum conservation law 
for the radiative decay in the wave also contains, along with the 4-momenta 
Pii P2-, Q, the 4-wavevector k. The process amplitude is calculated by standard 
Feynman rules, in which for the propagators of intermediate fermions (see Ap- 
pendix B) exact solutions are used of the corresponding wave equations. In this 
section we consider the external field of a monochromatic circularly polarized 
wave with the 4-potential 

= cos <p + £a 2/i sin ip, ip = kx, (40) 

where k fi = (ui,k) is the 4-wavevector, k 2 = (aik) = (a^fc) = (01^2) = 0, 
a\ = a\ = a 2 ; the parameter ^ = ±1 indicates the direction of the circular 
polarization (left- or rightward). Note that vectors a± : a<i and k form a right- 
handed coordinate system. The ^-matrix element of the given process, just as 
previously, can be represented in the following form: 

S = S + AS, (41) 

where So is the matrix element of the radiative decay of the massive neutrino 
in a vacuum and AS is the contribution induced by the wave field, 
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AS= /9F J [2 Sv nnV ^ A^ 4 )(nfc+ Pl -p 2 -g) (42) 
y/2EiV ■ 2E 2 V ■ 2q V n= _ 2 

Here pi, p 2 , q and E\, E 2l qo are the 4-momenta and the energies of the initial, 
final neutrinos and photon, respectively, n = 0, ±1, ±2 is the difference between 
the numbers of absorbed and emitted photons of the wave field. 

Note, that the matrix element of some process in the field of an electromagne- 
tic wave has usually the form of summation of n type (ffil), where — oo < n < oo 
2M . That only five values of n in our case are possible is extraordinary and is 



due to the following reasons. The process V{ — > is local with the typical scale 
Ax < 1/rrif (rrif is the mass of the virtual fermion). In this case the angular 
momentum conservation degenerates to spin conservation. Since the total spin 
of the particles participating in this process is no greater than 2, |n| max = 2 is 
the maximum difference between the numbers of absorbed and emitted photons 
of the external field (the photons of a monochromatic circularly polarized wave 
have a definite spin £ = ±1). The direct calculation supports this conclusion. 
A similar phenomenon has been discovered before [|I8j in studies of the effect 



of a circularly polarized wave on flavour - changing transitions of the massive 
neutrinos V\ <-> Vj (i ^ j) with |n| max = 1. 

Notice that in the uniform constant fields the decay V{ — > Vfl with rrii > rrij 
is valid only. This is due to the fact that the energy-momentum conservation 
law in this fields coincides with the one in vacuum. On the other hand, as it 
follows from Eq.(|E2f), the external electromagnetic wave field can also induce 



radiative transition with rrii < rrij forbidden without the field. Indeed, from 
the energy-momentum conservation law in the wave field 

nk + pi = P2 + q 

the relation follows 

m| — mj > —2n(kpi) 

In such a manner, the radiation decay V\ — ► with m a < rrij is possible on 
condition that n > 0. 

In the lowest order of perturbation theory, a matrix element of the radiative 
transition U{ — > Vj^ is described by the effective diagram represented in Fig. 3, 
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in which not only VF-boson, but also (for i = j) Z-boson is exchanged at the 
point x. The expressions for the fermion propagator in the wave's field, as well 
as some details of the awkward calculation of the invariant amplitudes J\4^ 
are given in Appendix B. Here we only present the final result: 



X 



M {a) = 



X 



+ 



G 



F 



16tt 2 (kq) 

{ J2(K u K* e + -Sijgi) ^\m\ [(jftyJ^mi) + (jfk)J 2 (m e ) 

£ 1 

3 _ 

^ijJ2Qq9q^ 2 q rn 2 q [(jfk)Ji(m q ) + (jfk)J 2 (m q )} }, 

Z q 

e G F ea(fF*) 1 



.2^2 



-%^(2T 3g Q 

l q 



(kq) 

((pi-P2)j)4 a \m q )+8 



me 



(43) 



m ^(jk)jt\m q ) 



(kq) 



e G F e 2 (fF°)(jF°q) 



16tt 2 



(kq)' 



1 



X { T,( K U K j£ + o^j^O J 5 ' ( '» ( ) 



e 2 
knQj), ^ kpo!). ^ , 



p/ = 2T 3/ - 4Q/ sin 2 



M (oi + i£(Ja 2 )^ 
f = t,Q. 



Here cr is the sign of the summation index n in Eq. fl42p (a = ±1), the index 
I indicates charged leptons (I = e, /i, r) and the index g indicates quark flavors 
(q = u, c, £, s, 6), T 3 y is the third component of the weak isospin and Qf is the 
electric charge in units of the elementary charge, mg and m q are the masses of 
the virtual leptons and quarks, Kij is the unitary lepton mixing matrix, 



Ji(rrif) = J dy J dr 



1 - W 



Jo 
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J. 



J 2 (m/) 
m/) 



= / dy 



o 



, 1 + y 2 

i-r 



00 







to 



0G 



% J dr(j 
(i o 



OT 



(44) 



J 



to, 



J. 



to 



m/) 



1 7 OO 



(77 



I OO 

% / drr(jQ+j'f)e 
o o 



z($(m/) 



2<xt) 



$(m/) = 



4r 



m 



1 



2/ 2 (M 



Jo, 



2 2 

2 fV 

2 ' 

mi- 



Here jo = sin r/r, ji = —djo/dr are so called Bessel spherical harmonics; the 
other denotations and quantities in (|43| ) and ([44]) have been introduced above. 
It is easy to see that the amplitudes MS n > are explicitly gauge-invariant and 
do not contain divergences. Note that the expressions we have obtained for 
the amplitudes A4^ of the radiative transition v\ —> Uj^ in the field of a 
monochromatic wave allow a simple check. Indeed, if in Eqs. ([43]) and ([44]) the 

— * — > 

wave frequency u tends to zero (i.e. /c M — > 0), the strengths £ and B of the 
electric and magnetic fields being fixed, then in this limit the amplitude of the 
decay V{ — ► Vfj in crossed field must be obtained (see Eqs. (P75|)— (P~5|) ) : 



AM[Eq. (13)] = M {n) 

n=-2 



k,. 



(0,0) 



(45) 



S, B - fix 

To prove the conclusion (|45|), it is necessary to take into account that: 
1) in the above limit, the tensor Fjfi) is expressed in terms of the strength tensor 



F^y and dual tensor F^ v of the crossed field: 



Fto 



k - 

S, B 



= i(j(Fn V + ie^F, 



fJLVJl 



(46) 



fix 



2) if we take advantage of the identity 
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(Aii 2 )a/3 + (A 2 il) 



1 



a/3 



a/3? 



(^4l^2)a/3 = Ai ap A2 p [3, (AiA 2 ) = ^4lap^2pa, 

where ^4.1^, are arbitrary antisymmetric 4-tensors, then the wave intensi- 
ty parameter sej (^Dj) can be related to the dynamic parameter xj (HO) by the 
expression: 



ml 



(47) 



3) to remove the undeterminacy arising at k — > in the expressions ( fffi ) for 
it is useful to take advantage of the limit relationship: 



k f 



(ka) 



(bFF) 



a 



(48) 



k^o (bFFaY 

where a and b are arbitrary 4- vectors (however, 6 is a timelike vector), 

4) in the crossed field limit the terms in (fJ3|) proportional to Sij do not 

make any contribution (because of kinematics reasons). 

In passing to the limit using fl35 Q -( |38|) the result ( ^5|) is reproduced immediately. 
The probability of the transition — > vfj in the wave field 



+2 

w = ^2 w 

n=-2 



(n) 



(49) 



is, in general, rather awkward. We shall give it only in the most interesting, 
from the physical point of view, case of the transition V{ —> Vjj with the initial 
neutrino Vi being ultrarelativistic [E v ^> m v ). Despite the fact that the wave 
intensity parameter aef (|39j) under laboratory conditions cannot be great (e.g. 
for laser fields uj ~ leV, £ < 10 9 V/cm, sel < 10~ 3 ), substantial enhancement 
of the transition probability is possible. The main effect of the enhancement 
connects with the decrease and even complete disappearance (for n > as we 
shall see below) of the suppression factor caused by the smallness of the neutrino 
mass. Recall that an analogous result was obtained for uniform and constant 
fields and was discussed in the above sections (see Eqs. flZ51)-(P7f) and (ffijp). In 
the ultrarelativistic limit, the probabilities at n < remain suppressed: 
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E v w K ' ~ U\a - A — ae 



mz 



e ' 



E^ -1 * ~ O ^ 



G 2 F ml 2 



(50) 



O (aG 2 F mlmlsef) . 



The other probabilities w^ (n = +1, +2) in the limit E v ^> m v are substantially 
simplified to be represented in the following form: 



7T 7T° J Z 



X 



-1 



dx 



"e e i 
1 — X 



1 (1 + x) 2 

»2 



Ui +1) (m 



(51) 



47T 7T° J Z 

'-(! + {) 



X 



fi dx 
i 1 + x 



'1-0 



+ 



'1 - .r ; 



|i +1 W 



Here, the integration is performed with respect to x = cost?, $ being the angle 
between the photon momentum <f and the wave vector k in the center-of-mass 
of the final neutrino vj and photon 7. Consequently, in the ultrarelativistic 
case in the integrals J5 (see Eqs. (@)) the substitution is needed 

(qk) ~ (1 + x)(p\k)/2. The comparison of the probabilities (|5U|)-(|5T|) of 
the transition vi — > ^7 in the wave field and the probability wo ( g0|) of the decay 
v% ~^ v il in vacuum shows that the probabilities w^ at n = +1, +2 do not 
contain suppression associated with the smallness of the neutrino mass (recall 
that the probability of the decay of an ultrarelativistic neutrino in vacuum is 
wo ~ m%jE v ). Below we estimate the ratio of the probability w (|9|) of the 
transition V{ — * Vj^ for a neutrino from a high-energy accelerator in the wave 
field of the laser type and the probability wq ( |2D|) of the decay in vacuum: 



wo 



10 



leVV E„u 



m v 



mi 



ioV x2 



(52) 
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where the wave intensity parameter ee 2 , for the laser type fields is: 



Such strong enhancement of the V{ —> vfj transition probability, even at rela- 
tively small wave intensity (ee 2 , — 10~ 3 ) appears rather impressive. 

The results obtained in this section may be of interest for astrophysics and 
cosmology. In particular, in the wave field the process of the photon splitting 
into the neutrino pair 7 — ► ViVj becomes possible. This process probability has 
the form: 



67T 07T° qo J 2 



+ (qk)\K ie K] e ^8 ij \ 2 \4 +1 \m e )\ 2 }. (54) 



As is easily seen from (|54j), this process probability, in the same way, is not 
suppressed by the smallness of the neutrino mass. It can be treated as an 
additional mechanism of the energy loss by stars. 

5 Conclusion 

In this work, in the framework of the Standard Model with fermion mixing, we 
have investigated the effect on the process z/j — ► vff of the massive neutrino 
radiative decay of external electromagnetic fields of various configurations: con- 
stant crossed field (section 2), constant uniform magnetic field (section 3), plane 
monochromatic wave's field (section 4). The analysis of the decay amplitudes 
and probabilities obtained leads to the following conclusion, which is the same 
for all the field configurations covered: an external electromagnetic field catal- 
yses the massive neutrino radiative decay. An especially strong enhancement 
occurs in the case of the ultrarelativistic neutrino radiative decay, since in this 
case the decay probability does not contain suppression caused by the smallness 
of the neutrino's mass. 

In section 2 (see eq. (|2g|)) we have estimated a number of gamma-quanta 



which could be resulted as the neutrino decay in a volume filled with a magnetic 
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field. Let us give here the estimation in the case of limiting in the laboratory 
conditions values of B and W (B ~ 10 9 G, W ~ 10 19 GeV): 



AW ~ 10- 2 ( JL;) 2 ( —) \K ie K 



10 9 G7 Vim, 



It is worth noting that the estimation we have presented is numerically small 
and seems likely that there is no possibility to carry out such neutrino experi- 
ment in the near future. 

Nevertheless, the mensioned above mechanism of the electromagnetic cataly- 
sis of the massive neutrino radiative decay is of interest in astrophysics where gi- 
gantic neutrino fluxes and strong magnetic fields can take place simultaneously 
(a process of a coalescence of neutron stars [p5fl , an explosion of a supernova 
of the type SN 1987A [j2§]). Let us estimate a relative flux of gamma-quanta 
which traverses a domain filled with a strong magnetic field (B ^> B e ): 

where Q v is the neutrino flux with the average energy E v ~ 1 MeV transversal 
to the magnetic field strength, L is the characteristic dimension of the domain. 
Gamma-quanta produced from the neutrino decay can be observed in astrophys- 
ical experiments provided that the domains with such strong magnetic fields 
exist. 

On the other hand, the results presented in the section 3 are, in our opinion, 
of interest in the cosmology of the early Universe. Indeed, in the recent papers 
it was pointed out a possibility of the generation of primordial strong magnetic 
fields through thermal fluctuation in the primordial plasma with the magnetic 
field strengths of order of 10 12 + 10 15 G J27| or 10 13 + 10 18 G |Z8j and coherence 
lengths of order of 10 -r- 100 cm. Let us estimate the neutrino lifetime in the 
case of existence of primordial small scale magnetic field strengths of order of 
~ 10 17 -T- 10 15 G. For this purpose we use the expression ( g7p we have obtained 
for the moderate energy neutrino (E v ~ kT ~ 1 MeV) lifetime in a strong 
magnetic field (B B e ) and get the following estimation: ~ 0.1 100 sec. 
This may be of interest in connection with the cosmological problem, conserning 
the contradiction between the COBE data on the cosmic microwave background 
anisotropics and the observed power spectrum of the large-scale structure Q29 
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A J a p calculation in constant 
electromagnetic field 

The amplitude corresponding to the diagram in Fig. 3 is calculated according to 
the conventional Feinman rules. In doing so, for the propagators of intermediate 
charged leptons exact solutions are used of the corresponding wave equations 
in the constant electromagnetic field. With the crossed field, the propagator of 
the charged lepton S^ F \x : y) in the proper time formalism |30] has the form: 



as 



S(X) = 



167T 2 J Q s- 



1 1 P ~ 

-(A 7 ) + -(AF 7 ) 75 
zs 2 




SC (XFF 7 ) + m-— ( 7 F 7 ) 



m 2 s + ^X 2 + ^-(XFFX) 
4s 12 v ' 



(55) 



(56) 



where X^ = 



x 



v) 



/i5 F [W , F^y are the field tensor and field dual tensor, e > is 



the elementary charge, 7/x , 7 5 are Dirac 7 -matrices (the metric, the convensional 



representation of Dirac 7 -matrices, etc. correspond to the book |23j| ), m is the 



mass of the charge lepton, the phase $(x, y) is determined in the following way: 
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X 



y 

KM) = MO + \f^-vV (57) 

Owing to d^K v — d v K^ = 0, the path of integration from y to x in ( |57j) is 
arbitrary and, therefore, 



®(x, y) + a;) = 0. (58) 

Using fl55p, the integration of J a/ 3 with respect to x and y (see Fig. 3) can 
be easily redused to an integration with respect to X = x — y (see Eq. (pT2|) ) . 
From (jI2f) and ( |55f) it is clear that the integrals with respect to X are Gaussian, 
so that they can be readily calculated: 



G = jd'Xe~ l ^ XRX + qX ) = -(4 7 r) 2 (deti?)- 1 / 2 e^^), 
G> = /^^ C -(^ M + ^)=tg, (59) 



d z G 



dq^dq v 



In the remaining double integral with respect to the proper times si, S2, it is 
convenient to pass to the dimensionless variables z, t: 



z = m 2 (si + S2), £ = , ds\ds2 = -rzdzdt, 

si + S2 2m 4 

< z < 00, — 1 < t < 1, (60) 

The substitution into the amplitude (|TT1 ) of the expression for J a p in the form 
of a double integral with respect to z, t results in the final expression 
and (ITJ). 

— > 

In the case of a constant uniform magnetic field B the propagator of the 
charged lepton S^ B \x, y) in the proper time formalism has the form: 
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S {B) (x,y) = 



i(3 



ds 



s sin(/?s) 



x 



cos((3s)(XAj) + zsin(^s)(X^7)7 5 



sin(/3s) 



(XA7) + m [2 cos(/3s) - sin(^s)(7^7)] 



(61) 



x exp 



m s + 



4s 4tan(/3s 



-(XAX) 



where ip^ = F^/B, Cp iW = F^/B are the dimensionless field tensor and dual 
field tensor of the constant magnetic field, A a p = ((p(p) a p, A a p = (<p0) a p, 
(3 = eB, X^ = (x — the phase $(x,y) is described in (|57|). Note that 



the propagator (£JT]) can be represented in a fully covariant form, because the 



parameter (3 in a purely magnetic field can be rewritten as (3 = eB = \j—F 2 /2. 
The calculation procedure for the tensor J a p ( \T2\ ) in the case of a constant mag- 
netic field, though more awkward, does not, in principle, differ from the case of 
a crossed field. 



B On the calculation of the /S-matrix element in the field 
of a monochromatic circularly 
polarized wave 

The propagator of the charged fermion in the field of a plane wave with a 4- 
potential = A^ip) of the general form can be obtained by method, given in 
|3Tj| , and has the form: 



S(x,y) = 



A / \ / 1 

d p ( ejkA \ p + ruf ( ejAk 
(2tt) 4 [ ~ 2(kp) ) p 2 -m)\ ~ 2{kp) 



(62) 



x exp |i - p(y -x) + j^-j J dip {e f {pA) + ^e 2 f A 2 ) } 
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where = A^(p), p = kx, A'^ = A^(p'), p' = ky, k is the 4-wavevector 
(k 2 = 0), ej and mf are the charge and the mass of fermion, respectively. 
In the case of the circularly polarized wave with 4-potential 

A fl (p) = a lfX cos(p) + a 2M sin(^), 

where 4- vectors and are orthogonal to the 4-wavevector k^. 

(0,10,2) = (a\k) = (02k) = 

the expression for propagator may be represented in the following form: 



(63) 



(64) 



S(x,y) = 



d 4 p 



1 



efkA'\ p + rrif 



x 



» 4 
exp ji — p(y — x 
e f (a 2 p) 



2(kp) J p 2 — m 2 
e f (a lP ) 



1 



ejAk 



(kp) 



2(kp) 
(sin(<//) — sin((/?)) 



(65) 



2 2 

Since the power of the exponent of the propagator ( |55| ) contains nonlinear 
functions of coordinates x and y (sin(<p), sin(<//), cos(^) and cos(p')), it is 
convenient to expand the corresponding part of an exponent in the Fourier 
expansion with the coefficients of expansion been proportional to the Bessel 
functions |2Jj . 

Given the integration over one of the momenta in the loop d A q it is convenient 
to use the following basis: p M , h^p^ h^ v p v , h^h^pp, where 



h^v — k^aiv k v o\^, h^y — k^02 v k^02^,, 
By using the known relation 



(66) 



Jo (Vb 2 -2bccosa + c 2 ) = J (b)J (c) + 2 £ J s (b)J s (c)cossa, (67) 
^ ' s=l 

one can reduce the remaining infinite series to five terms, which may be brought 
to the form ([42]). 
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